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ABSTRACT 


There  are  some  indications  that  in  large  measure  the  density  anomaly  field  of  the 
California  Current  System  (CCS)  follows  temperature,  salinity  playing  a  significantly 
lesser  role.  These  indications  have  been  used  as  justification  for  ignoring  salinity  vari¬ 
ations  in  dynamic  models.  An  extensive  data  base  of  simultaneous  temperature-salinity 
observations  taken  in  the  CCS  is  used  to  calculate  and  assess  the  saline  contribution  to 
specific  volume  anomaly,  dynamic  height,  and  density  variance.  The  results  show  that 
the  distribution  of  salinity  can  be  important  in  defining  the  large  scale  circulation  of  the 
CCS,  and  that  the  local  variability  can  be  quite  significant. 
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I.  INTRODUCTION:  IMPORTANCE  OF  SALINITY 


At  a  given  pressure,  the  density  of  seawater  depends  on  temperature  and  salinity. 
Except  for  high  latitude  regions  and  near  sources  and  sinks  of  salt,  teir  erature  is  usually 
considered  to  vary  more  than  salinity  in  the  upper  ocean,  and  is  generally  considered  to 
be  the  major  source  of  density  variations  (Pickard  and  Emery,  1982).  Recently  Cooper 
(1988)  has  demonstrated  the  importance  of  salinity  as  a  consideration  in  modeling 
motions  in  the  tropics.  This  result  is  surprising  because  in  the  tropics  the  thermal  ex¬ 
pansion  coefficient,  which  increases  with  temperature  and  salinity,  is  large,  and  the 
salinity  contraction  coefficient,  which  decreases  as  temperature  and  salinity  increase,  is 
relatively  small,  and  temperature  variation  might  reasonably  be  expected  to  dominate 
forcing  of  density  differentials.  Is  it  possible  that  salinity  variability  in  the  California 
Current  System  (CCS),  which  is  geographically  between  the  polar  regions  and  the  trop¬ 
ics,  is  also  an  important  consideration  for  modeling  dynamics?  A  close  examination  of 
the  southern  CCS  area  in  Figure  I  suggests  that  salinity  variation  could  have  as  im¬ 
portant  an  effect  in  the  CCS  as  in  the  tropics. 

The  CCS  is  in  the  eastern  quadrant  of  the  North  Pacific  anticyclonic  gyre,  and,  ac¬ 
cording  to  Hickey  (1979),  is  composed  of  four  currents.  The  California  Current  (CC)  is 
on  the  surface  and  extends  as  far  as  1000  km  from  the  coast.  A  surface  poleward  flow 
exists  during  the  fall  and  winter  north  of  Point  Conception  called  the  Davidson  Current. 
South  of  Point  Conception,  and  inshore  of  the  Channel  Islands  is  another  surface  flow 
called  the  Southern  California  Countercurrent.  Finally,  the  California  Undercurrent 
(CU)  is  a  subsurface  poleward  flow  that  occurs  over  the  continental  slope.  Typical  mean 
velocities  associated  with  the  CC  are  on  the  order  of  10  cm  s-',  but  superimposed  on  the 
mean  flow  are  small  scale  eddies  and  jets  (Bernstein  et  al,  1977;  Mooers  and  Robinson, 
1984;  Rienecker  et  al.,  1985,  1988).  The  meanders  can  have  wavelengths  of  several 
hundred  kilometers,  can  intensify  and  separate  from  the  mean  flow  and  can  exist  as 
isolated  cyclonic  or  anticyclonic  eddies  (Bernstein  et  al.,  1977).  Baroclinic  jets,  typically 
~70  km  wide,  with  peak  velocities  near  80  cm  s“‘  are  embedded  in  the  eddy  field  (Mooers 
and  Robinson,  1984;  Rienecker  et  al.,  1985;  Flament  et  al.,  1985;  Kosro  and  Huyer, 
1986).  Mean  surface  temperatures  in  the  CCS  range  roughly  from  10  to  20®  C  while 
surface  salinities  increase  from  32  psu  in  the  north  to  34  psu  in  the  south. 
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Figure  1.  Dynamic  Height  of  I 
Lcvitus  (1982)  data,  (f 
as  in  Lcvitus  data,  and 
tervals  arc  10  ifyn  an  ir 
1988). 


lie  Tropical  Ocean  0/400  in:  Calculated  from 
)  Salinity  set  to  constant  35,  (b)  salinity  varying 
(c)  diirorcncc  between  (a)  and  (b).  Contour  in- 
(a)  and  (b),  and  5  dynan  in  (c)  (From  Cooper, 
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A  typical  section  across  the  CCS,  from  San  Francisco  Bay  southwest  to  33®  N, 
134®  W  (Figure  2),  shows  that  isopycnals  very  closely  follow  isotherms,  and  intersect 
isohalines.  The  implication  is  that  temperature  variation  is  much  more  important  than 
salinity  variation  in  driving  the  density  gradient.  On  the  other  hand,  examination  of 
Figure  3,  which  depicts  the  coastal  area  of  the  CCS  just  north  of  Point  Arena  (at 
~39®  N  and  near  the  shoreward  end  of  Figure  2)  shows  that  at  a  depth  of  100  m,  about 
a  third  of  the  cross-shore  density  gradient  over  a  degree  of  longitude  is  due  to  salinity 
change.  In  the  CCS,  unlike  all  other  eastern  boundary  currents,  salinity  increases  with 
depth  (Wooster  and  Reid,  1963).  Therefore,  if  this  cold,  salty  100  m  water  were  upwelled 
to  the  surface,  the  resulting, horizontal  density  gradient  would  be  one  third  greater  than 
if  the  salinity  profile  were  constant. 

The  approximate  equation  of  state  for  seawater  is: 

p  =  Po[I-«(7'-ro)  +  /?(S-5o)],  (I) 

which  can  be  rewritten  as: 

Ap  =  -  (xpoAr  +  /IpoAS,  (2) 

where  a  and  P  are  the  expansion  and  contraction  coefTicients  for  temperature,  T,  and 
salinity,  S,  respectively,  p  is  density,  S®  is  the  reference  salinity  (35  psu)  and  To  is  the 
reference  temperature  (0®  C). 

This  equation  can  be  employed  to  get  a  rough  idea  of  the  large  scale  importance 
of  salinity  if  typical  CCS  values  for  T  and  S  are  used.  Entering  the  International 
Oceanographic  Tables  (UNESCO,  1987)  for  a  and  p  with  15®C  and  33  psu,  and  con¬ 
sidering  typical  surface  temperature  and  salinity  ranges  of  10®  C  and  2  psu  for  the  CCS, 
computations  of  the  contributions  of  temperature  and  salinity  to  the  density  anomaly 
are  10.5  and  7.5  kgm-\  respectively.  Annually,  a  mean  temperature  of  12®  C  plus  or 
minus  3®  C  variation,  associated  with  a  mean  salinity  of  33  with  .3  psu  variation,  is 
roughly  characteristic  of  the  central  region  of  the  CCS.  Use  of  the  1987  tables  here 
shows  a  1.7  and  2.7  kgtrr^  contribution  to  the  annual  density  anomaly  change  for 
salinity  and  temperature,  respectively,  over  the  course  of  a  year.  Vertically  a  mean 
temperature  at  500  m  of  about  5®  C  plus  or  minus  5®  C  surface  to  1000  m,  in  conjunction 
with  an  appro.ximate  mean  salinity  of  34,  plus  or  minus  0.5  psu,  yields  density  anomaly 
changes  of  7.2  kg  nr^  from  temperature  and  3.8  kg  nr^  from  salinity. 
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These  crude  calculations  seem  to  indicate  that  the  importance  of  salinity  variation 
should  not  be  summarily  dismissed  when  considering  CCS  density  gradients  in  any  di¬ 
mension.  Yet  often  dynamic  models  of  the  CCS  (e.g.,  Batteen  et  a/., 1989)  do  not  allow 
salini,t}'^contributions  to  density  gradients.  What  is  the  price  of  this  omission? 

The  objective  of  this  thesis  is  to  address  that  question  by  investigating  the  saline 
impact  on  the  CCS.  The  data  set  used  is  described  in  Section  II.  A  general  description 
of  the  characteristic  salinity  of  the  CCS  is  provided  in  Section  III.  In  Section  IV  several 
techniques  are  used  to  analyze  the  data.  Mean  values  for  each  of  the  constituents  of  the 
thermosteric  anomaly,  5,d,  and  5,„  are  computed  at  various  depths  and  analyzed.  Mean 
dynamic  heights  are  computed  from  salinity  and  temperature  observations  and  are 
compared  to  dynamic  heights  computed  from  temperature  and  constant  salinity  values. 
The  equation  of  state  for  seawater  is  differentiated  with  respect  to  temperature  and 
salinity,  variances  and  covariances  of  temperature  and  salinity  are  computed,  and  the 
contributions  of  each  to  variance  in  the  density  field  are  determined.  A  summary  of  the 
results  is  presented  in  Section  V  and  recommendations  for  future  work  are  in  Section 
VI. 
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Figure  2.  CalCOFI  Line  60,  July:  33*  N,  134* W,  to  37*  N,  123*  W:  (a)  I  cni- 
pcraturc  (*C),  (b)  Salinity  (ppt),  (c)  Density  Anomaly  (from  Lynn  et  at.. 
1982). 
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Figure  3.  Temperature,  Salinity,  and  Density  Gradients:  Near  Pt.  Arena  (39°  N), 
at  100  ni  depth:  (a)  Temperature  (°  C),  (b)  Salinity  (ppl),  (c)  Density 
Anomaly  (from  lluyer  and  Kosro,  1987). 
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II.  DATA  DESCRIPTION 


A.  SOURCES  AND  TYPES  OF  DATA 

The  data  base  used  for  all  calculations  was  compiled  by  Dynalysis  of  Princeton,  and 
is  described  in  detail  in  their  final  report  to  the  U.S.  Department  of  the  Interior 
(Blumberg  et  ah,  1984).  Coverage  extends  along  the  west  coast  of  North  America  to 
133  °  W,  and  from  20®  to  50®  X.  Sources  include  the  National  Oceanographic  Data 
Center  (NODC)  historical  data  archives,  the  Fleet  Numerical  Oceanography  Center 
(FNOC)  data  files,  California  Cooperative  Oceanic  Fisheries  Investigation  (CalCOFI) 
surveys  which  have  been  made  systematically  since  1949,  and  the  Coastal  Ocean  Dy¬ 
namics  Experiment  (CODE).  CODE  data  are  the  most  current,  dating  from  as  recently 
as  1983,  and  cover  the  coastal  area  near  Point  Arena  within  a  one  degree  square  domain. 
The  highest  density  of  NODC  data  is  the  result  of  observations  from  Oregon  State 
University  standard  surveys  along  the  "Newport"  line  (~  45®  N)  and  University  of 
Washington  surveys  of  the  Columbia  River  plume  (46-47®  N).  Original  files  overlapped, 
but  the  data  have  been  carefully  screened  by  Dynalysis  to  eliminate  duplicate  observa¬ 
tions.  The  combined  data  set  consists  of  over  300,000  data  points  with  observations  at 
up  to  33  standard  levels  to  4000  m.  The  more  than  40,000  of  these  points  which  included 
any  simultaneous  temperature  and  salinity  observations  at  500  m  or  above,  were  con¬ 
sidered  here,  w’hile  those  including  only  temperature  were  not  used.  Original  observa¬ 
tions  of  salinity  were  reported  in  units  of  parts  per  thousand  (ppt)  rather  than  practical 
salinity  units  {psu).  However,  because  the  difference  in  units  is  very  slight  (Lewds  and 
Perkin,  1981)  and  not  significant  at  the  scales  discussed  here,  the  currently  preferred  psu 
is  used  throughout  this  thesis.  Most  of  the  retained  records  are  from  CalCOFI  surveys, 
which  were  collected  at  average  intervals  of  74  km  along  lines  spaced  74  km  apart  cov¬ 
ering  from  San  Francisco  to  Baja  California  (Reid,  1988).  As  a  result  the  CalCOFI 
collection  defines  for  the  most  part  the  maximum  resolution  possible  from  the  data  set, 
exceptions  being  along  the  Oregon  coast,  and  near  Point  Arena,  w’here  finer  resolution 
may  be  possible.  Generally  the  number  of  observations  increases  toward  the  coast,  with 
greatest  concentration  along  the  coasts  of  California  and  Oregon.  There  are  very  few 
data  more  than  500  km  off  the  coast  or  between  38®  and  42®  N.  Discussion  here  is 
limited  to  the  area  bounded  by  the  extent  of  this  data. 
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B,  DATA  MANIPULATION 

For  this  study,  data  were  placed  in  bins  identified  by  one  degree  squares  centered 
at  half  degree  latitude-longitude  intersections,  by  month  or  season,  and  by  depth  level. 
Information  so  cataloged  is  referred  to  as  a  data  block  throughout  this  paper.  Standard 
depth  levels  0,  10,  20,  30,  50,  75,  100,  125,  150,  200,  250,  300,  400,  and  500  m  were 
considered.  Choice  of  500  m  as  the  maximum  depth  level  was  made  due  to  data  avail¬ 
ability  and  because  the  phenomenon  of  interest,  the  CCS,  has  been  classically  considered 
shoaler  than  that  depth,  at  least  within  300  km  of  shore  (Hickey,  1979).  Choice  of  one 
degree  squares  was  made  again  due  to  data  availability;  the  CalCOFI  sampling  grid  does 
not  support  finer  resolution. 

Plots  of  the  numbers  of  observations  associated  with  each  data  block  for  the  months 
January,  April,  July  and  October  are  in  Appendix  A.  When  mean  values  and  variances 
were  calculated,  at  least  ten  observations  in  a  data  block  were  required,  except  in  the 
calculation  of  mean  dynamic  height,  where  five  observations  were  used.  Mean  values 
generated  from  these  numbers  of  observations  were  considered  statistically  significant 
and  used  unweighted  in  subsequent  analyses.  Weighting  based  on  redundancy  was  not 
applied  to  avoid  biasing  results  in  favor  of  areas  of  frequent  observation.  The  determi¬ 
nation  of  ten  observations  as  statistically  significant  was  done  subjectively  after  studying 
the  data.  Unfortunately  specific  years  over  which  statistics  for  a  given  data  block  were 
not  considered,  and  a  bias  toward  a  particular  year,  or  years,  may  well  exist.  For  mean 
dynamic  height  calculations,  only  casts  which  include  temperature  and  salinity  data  at 
all  standard  depths  from  the  surface  to  500  m  were  used. 
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III.  CHARACTERISTIC  SALINITY  OF  THE  CALIFORNIA  CURRENT 

SYSTEM 


A.  GENERAL  FEATURES 

Along  .the  west  coast  of  the  United  States,  sea  surface  salinity  generally  increases 
from  about  32  psu  near  the  Washington-Canadian  border  to  about  34  psu  near  Baja 
California.  At  a  given  latitude  salinity  tends  to  be  maximum  at  the  coast,  decrease  off¬ 
shore,  then  increase  again  at  the  western  limits  of  the  CCS  (Figure  4).  That  is,  the 
isohalines  describe  a  trough,  the  axis  of  which  extends  from  the  coast  near  the  mouth 
of  the  Columbia  River  near  46°  15'  N  to  a  point  about  500  km  off  the  coast  at  Baja 
California.  In  the  southwest  corner  of  the  CCS,  approximately  700  km  off  the  coast,  the 
isohalines  are  oriented  northwest-southeast  with  a  relatively  tight  gradient,  of  the  order 
0.1  psu  in  40  km. 

As  mentioned  earlier,  the  stratification  of  the  CCS,  unlike  that  of  other  eastern 
boundary  currents,  has  salinity  generally  increasing  with  depth.  Profiles  (Figure  5)  are 
marked  by  permanent  haloclines:  sharp  and  shallow,  150  m  (Huyer,  1983)  in  the  north 
and  less  distinct  and  deeper,  200  m,  to  the  south  (Emery  and  Dewar,  1982).  .At  500  m 
the  mean  salinity  increases  from  34.1  psu  near  50°  N  to  34.5  psu  near  20°  N. 

B.  VARIABILITY 

The  most  dramatic  seasonal  change  to  the  surface  salinity  in  the  CCS  results  from 
river  discharge.  During  the  summer  the  only  important  sources  of  fresh  water  runoff 
along  the  West  Coast  are  the  Fraser  River,  which  discharges  through  the  Straits  of 
Georgia  and  Juan  de  Fuca  at  48°  30  ’  N,  and  the  Columbia  River  at  46°  15'  N.  The 
Fraser's  fresh  water  is  mixed  deep  and  its  signature  lost,  but  the  direct  effect  of  the 
Columbia  is  visible  as  far  south  as  San  Francisco  (Huyer,  1983).  During  spring  an  area 
of  tight  isohaline  spacing,  roughly  0.1  psu  in  30  km,  exists  within  ~500  km  radius  of  the 
Columbia's  mouth.  Through  the  course  of  a  year  the  area  of  tight  gradient  (Figure  4) 
gradually  decreases  until  by  winter  the  strong  gradient  exists  only  in  the  immediate  vi¬ 
cinity  of  the  mouth  of  the  river  and  northward  along  the  coast  of  Washington. 

Using  CalCOFI  data  collected  as  far  north  as  San  Francisco  through  1978,  Lynn 
and  Simpson  (1987)  have  e,\amined  the  physical  characteristics  of  the  CCS  including 
salinity.  The  salinity  minimum  is  at  the  surface  off  the  coast  of  northern  California  and 
50  to  100  m  deep  in  the  Baja  California  region.  Near  southern  California  there  is 
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arc  ppt.  (a)  Sununcr  (b)  Winter 


generally  a  minimum  at  the  surface,  but  in  the  autumn  a  subsurface  minimum  appears. 
Lynn  and  Simpson  (1987)  have  noted  significant  correlations  between  the  location  of  the 
salinity  minimum  and  the  behavior  of  the  CC.  In  January,  a  strong  inshore  poleward 
current  exists  while  the  salinity  minimum  near  Baja  California  is  below  the  surface  and 
offshore  of  the  continental  slope.  As  this  minimum  moves  inshore  and  to  the  surface, 
the  inshore  flow  becomes  strongly  equatorward.  In  October,  the  salinity  distribution  is 
similar  to  that  of  January,  but  the  inshore  flow  is  equatorward.  Elsewhere  in  the 
CalCOFI  region,  the  area  of  minimum  salinity  is  consistently  near,  but  seaward  of,  the 
maximum  equatorward  flow  of  4  to  20  cm  s-'.  Because  of  the  high  correlation  of  salinity 
minimum  to  velocity  maximum  to  a  depth  of  100  m,  Lynn  and  Simpson  (1987)  suggest 
that  the  low  salinity  core  can  describe  the  path  of  the  CC. 
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IV.  EFFECT  OF  SALINITY 


A.  STRATEGY 

Several  techniques  were  used  to  analyze  the  data  and  are  described  in  the  following 
sections:  specific  volume  anomaly  (Section  I VB),  variance  and  correlation  of  the  specific 
volume  anomaly  constituents  (Section  I  VC),  mean  differences  of  the  specific  volume 
anomaly  constituents  (Section  IVD),  dynamic  height  (Section  IVE),  and  propagation 
of  the  temperature-salinity  variance  (Section  IVF).  Each  method  was  applied  at  depths 
10,  100,  200,  and  500  m.  Except  in  Section  IVF,  twelve  mean  months  were  examined. 
In  Section  IVF  the  statistical  techniques  invoked  called  for  use  of  maximum  number  of 
observations,  so  the  data  were  partitioned  into  three-month  seasons;  winter,  for  exam¬ 
ple,  was  defined  as  January-March.  Graphics,  except  where  otherwise  noted,  were  gen¬ 
erated  for  all  four  depths,  each  time  period,  and  each  technique.  Practical  limits  on 
space  prevent  inclusion  of  all  these  plots  in  this  thesis.  l  Only  representative  examples, 
which  when  contonrod  were  hand  contoured,  are  reproduced  in  the  body  of  the  te.\t.  In 
particular,  the  month  of  July  is  often  chosen  as  the  illustrative  period  because  the  coastal 
upwelling  phenomenon  is  usually  at  its  peak  during  that  month. 

B.  SPECIFIC  VOLUME  ANOMALY 

1.  Method 

Specific  volume  anomaly  (defined:  V—  where  V  is  specific  volume,  and 
subscripts  in  the  second  term  indicate  salinity  35  psu,  temperature  0®  C,  and  pressure  in 
db)  can  be  expressed  as  a  sum  of  constituents  related  to  the  parameters  that  define  the 
density  of  sea  water: 


5  —  {5(  -f-  S/p  -f-  djp  Sjip, 

(1)  (2)  (3)  (4)  (5)  (6) 

Subscripts  t,  s,  and  p  stand  for  temperature,  salinicv  "nd  pressure,  respectively.  In  this 
linear  construct,  specific  volume  anomaly,  also  tc»T  ted  steric  anomaly,  is  the  sum  of  a 
part  due  to  variation  in  temperature  from  0*  C,  salinity  other  than  35  psu,  their  inter¬ 
action,  and  the  interaction  of  deviant  temperature  and  salinity  with  pressure  (Sverdrup 

I  Unpublished  data  are  on  file  at  the  Department  of  Oceanography,  Naval  Postgraduate 
School,  Monterey,  CA. 
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et  ai,  1942:57-59).  Since  the  term  <5„,  is  much  smaller  than  the  other  terms,  it  is  virtually 
always  ignored,  and  is  not  considered  in  this  thesis.  The  sum  of  terms  in  parenthesis 
represents  the  total  combined  contribution  of  temperature  and  salinity  calculated  at 
surface  pressure,  is  called  the  thermosteric  anomaly,  A,„  and  is  by  far  the  greatest  con¬ 
tribution  to  the  total  specific  volume  anomaly.  Sverdrup  (1933)  introduced  the  use  of 
the  thermosteric  anomaly,  which  is  easily  determined  in  total  from  the  measured  or 
computed  surface  density,  as  a  computational  convenience.  After  the  density  anomaly 
calculated  at  zero  pressure  is  converted  to  only  the  two  values  5,,,  and  5,,  are  neces¬ 
sary  to  determine  specific  volume  anomaly. 

Traditionally  then,  since  Sverdrup  (1933),  only  three  tables  (or  calculations)  are 
used  to  determine  S.  However,  it  can  sometimes  be  quite  useful  to  know  the  individual 
contributions  to  A„.  For  e.Kample,  if  5,  is  small  or  constant  while  5,  is  large  and  fluctu¬ 
ating,  a  dynamic  modeler  might  simplify  his  model  by  keeping  salinity  constant.  Simi¬ 
larly,  if  a  feature  can  be  defined  by  its  salinity  signature,  an  observationalist  .Jght 
identify  the  feature's  boundary  by  a  sharp  gradient  in  5, . 

The  equation  of  state  for  seawater: 


where  k  is  the  secant  bulk  modulus  and  P  pressure,  translates  to  an  equation  for 
thermosteric  anomaly: 


-  9.7266204  x  10"^  (5) 

or 

A„- - ; - ; — ! - 5 -  9.7266204x  10“',  (6) 

(<>w  -  %/)  + 

/l*0 

where  is  the  density  of  reference  pure  water,  a,  b,  and  c  are  coefficients  in  power  se¬ 
ries  on  temperature,  and  is  a  coefficient  on  salinity  (UNESCO,  1987).  This  inverted 
sum  of  terms  in  S  and  T,  cannot  readily  be  analytically  sorted  into  parts  due  solely  to 
temperature,  salinity,  and  their  interaction.  Therefore,  some  other  means  of  defining 
d„  5„  and  8„  is  needed.  In  UNESCO  (1987)  equation  (5)  is  modeled  by: 
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(7) 


^stp—  ^stO  +  ~  ^^35/0)  +  (^jOp  ~ ‘^jOo) 

(I)  (2)  (3)  (4)  (5) 

{a)  (b)  (c) 

(Here  d  is  used  in  lieu  of  S  to  alleviate  confusion  that  might  otherwise  result  when 
comparing  terms  of  equation  (3)  to  terms  of  equation  (7)).  Clearly  specific  volume  cal¬ 
culated  for  a  given  salinity  and  temperature  at  surface  pressure,  d,^,  is  precisely  the 
thermosteric  anomaly.  Term  (b),  the  difference  between  d  calculated  at  temperature  and 
depth  with  salinity  35  psu  and  d  calculated  for  the  same  temperature  but  at  surface 
pressure,  describes  the  interplay  of  temperature  and  pressure,  <5,,  .  Similarly,  term  (c) 
represents  5,^.  (Terms  a,  b,  and  c  correspond  to  UNESCO,  1987,  tables  III,  V,  and  VI.) 

Now,  consistent  with  the  currently  accepted  equation  of  state  for  seawater 
(UNESCO,  1981),  how  does  one  model  the  individual  constituents  of  the  thermosteric 
anomaly?  Term  (3)  in  equation  (7)  is  the  value  of  S  that  would  result  at  zero  pressure 
if  salinity  were  a  constant  35  psu.  Therefore  {d,^  —  d^^)  represents  the  part  of  3  due  to 
salinity  deviation  from  35  psu  and  interaction  of  temperature  with  that  deviation  at  a 
certain  pressure  p  ,  or: 


^stp  -  <(35/0  +  ^ip  +  ^sp)' 


Similarly: 


^sip  ~  ^sOO  “  i^st  +  ^tp  +  ^sp  +  ^/)- 


Subtraction  of  equation  (9)  from  equation  (8)  yields: 


(8) 

(9) 


=»  d^sio  —  d^.  (10) 

It  would  appear  d^a  and  d^  can  serve  as  models  for  3,  and  3,.  Using  that  assumption, 
equation  (3)  can  be  reconfigured  from  appropriate  elements  of  equation  (7): 
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-  ^35(0 
=  ^jOO 


^st  -  ^st  —  ^4/0  ~  ^35t0  ~  ^sOO 

* 

^tp  =  ‘^iUp  ~  ^35t0 

(11) 

^sp  —  ^sOp  ~  ^sOO 

^stp  =  ^35iO  +  ^jOO  +  i^stO  ~  ^35ri)  ~  ^iOo)  +  (^35fp  “  ^35rt)  +  (^4op  “  ^^sOo)’ 


(1)  (2)  (3)  (4)  (5) 

Collection  of  terms  shows  equation  (1 1)  to  be  identical  to  equation  (7)  and  therefore  this 
model  to  be  consistent  with  the  UNESCO  (1987)  table  formulation.  Descriptively  this 
formulation  should  also  be  acceptable:  the  temperature  contribution  is  defined  as  the 
specific  volume  that  would  result  from  the  observed  temperature,  if  salinity  w’ere  35  psu 
and  pressure  0  db;  the  salinity  contribution  at  the  same  depth  is  similarly  the  specific 
volume  anomaly  that  would  result  from  actual  salinity  and  0®  C;  the  combined  term  is 
simply  the  difference  of  the  sum  of  those  two  terms  and  A„  or  the  amount  by  which  their 
sum  misrepresents  A„. 

It  follows  that  (5,  +  5,^)  can  be  represented  by  and  (5,  +  5  J  by  It  w'as 
decided  to  use  this  combined  temperature  and  temperature-pressure  term  in  order  to 
compare  total  saline  effect  to  temperature  effect  at  a  given  pressure.  As  a  notational 
convenience,  throughout  the  remainder  of  this  thesis  “3/  and  "3^  "  will  be  used  in  lieu 
of  "(5,  +  3,^) "  and  "{3,  +  3,,)  '  and  be  considered  synonymous  w'ith  "  and  In 
this  formulation,  the  term  5„  =  —  3r—  3s  will  include  pressure  effects.  Since  units 

used  to  describe  specific  volume  anomaly  will  in  every  case  be  10**  kg~\  they  w'ill 
hereafter  be  deleted  in  the  text. 

2.  Results 

For  each  month  and  the  depths  10,  100,  200,  and  500  m,  each  of  the  fields  3t, 
ds,  and  5„  were  computed  and  plotted  (see  note  1).  Generally,  the  specific  volume 
anomaly  field  consists  of  a  dome  centered  near  45“  N,  124“  W,  a  depression  centered 
near  San  Francisco,  and  a  ridge  parallel  and  adjacent  to  the  Baja  California  peninsula 
with  axis  south  to  north  about  3“  longitude  from  shore.  In  July  (Figure  6)  all  three 
features  are  w'ell  established.  The  Columbia  River  outflow  is  at  its  peak  during  late  June 
and  early  July  (Huyer,  1983),  and  the  warm,  fresh  surface  plume  results  in  enhancement 
of  the  northern  dome.  July  is  also  a  period  of  strong  coastal  upwelling  (Huyer,  1983), 
and  the  cold,  salty  upwelled  water  results  in  intensification  of  the  central  depression. 
The  horizontal  gradient  in  3  decreases  with  depth,  particularly  in  the  northern  region. 
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At  100  m  (Figure  7)  for  example,  there  is  no  evidence  of  the  warm,  fresh  core  so  pre¬ 
valent  on  the  surface.  At  500  m  (Figure  8)  the  gradient  is  very  weak  everywhere  with 
dense  cores  centered  about  2*  off  the  coast  at  San  Francisco  and  the  northern  Baja 
California  peninsula.  An  important  observation  is  that  d„  (Figure  9)  values  are  one  to 
two  orders  of  magnitude  less  than  the  sum  of  5r  (Figure  10)  and  ds  (Figure  11).  This 
shows  the  steric  anomaly  to  be  very  nearly  equal  to  the  sum  of  a  temperature  and 
salinity  part  so  that  this  formulation  is  an  effective  linearization.  In  all  cases  6„  is  neg¬ 
ative,  increasing  slightly  in  absolute  value  with  depth;  values  of  -3  to  -6  being  common 
at  10  m  and  -6  to  -9  more  typical  at  500  m  (Figure  12).  Gradients  in  are  everv'where 
small;  however  the  cross  term  becomes  more  negative  in  areas  where  temperature  and 
pressure  act  together  to  increase  S,  and  less  negative  when  temperature  and  salinity  act 
in  conjunction  to  decrease  S.  When  salinity  and  temperature  oppose  change  in  specific 
volume  anomaly,  5„  remains  nearly  constant;  the  temperature-pressure  interaction  acts 
to  mitigate  change  in  specific  volume  anomaly  that  results  from  changes  in  either. 

Advection  of  subarctic  water  south  from  the  point  where  the  North  Pacific 
Current  feeds  the  CC,  results  in  cold  (low  S,)  and  fresh  (high  S,)  tongue-like  structures 
superimposed  on  the  other  characteristic  features  of  the  CCS  described  in  part  below. 
In  the  immediate  vicinity  of  the  coast,  salinity  and  temperature  act  together  to  decrease 
steric  anomaly.  Near  the  surface  this  is  most  apparent  just  north  of  San  Francisco  Bay 
in  July  where  large  scale  upwelling  results  in  cold,  salty  water.  The  effects  on  density 
of  temperature  and  salinity  are  positively  correlated  with  specific  volume  anomaly  near 
45“  N,  127"  W,  where  the  warm,  fresh  discharge  of  the  Columbia  River  results  in  the 
geopotential  dome  described  above.  This  feature  is  strongest  in  sumj.ier,  weakest  in 
winter,  and  most  apparent  in  the  ds  field.  Such  seasonality  is  in  keeping  with  the 
Columbia  River's  outflow  cycle.  There  also  persists  a  fresh  dome  at  the  mouth  of  the 
Strait  of  Juan  de  Fuca,  but  it  is  restricted  to  the  immediate  vicinity  of  the  coast.  (Since 
the  fine  scale  salinity-temperature  interplay  associated  with  the  drainage  from 
Vancouver  Island  is  not  resolvable  on  the  coarse  scale  used  here,  activity  in  that  area  is 
generally  not  included  in  discussion  here.) 

Near  the  surface,  throughout  the  year,  mean  gradients  of  across  the  CCS 
(neglecting  small  scale  features  like  those  described  above)  are  at  least  50  percent  greater 
than  those  of  At  a  depth  of  10  m,  absolute  values  of  which  range  from  about  80 
to  400  north  to  south  through  the  year,  fluctuate  more  season  to  season  and  are  gener¬ 
ally  larger  than  values  of  dj  which  range  typically  from  50  to  250.  However  the  gradient 
in  5r  diminishes  more  rapidly  with  depth  than  does  that  of  5s  ;  the  result  is  that  at 
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Figure  7.  Steric  Aiiomnly,  July,  100  iii:  Units  arc  I0-* Ag  ',  and  tlic  contour 
interval  is  10  units. 
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Figure  8.  Steric  Anomaly,  July,  500  m:  Units  arc  10  »/n’Ag-',  and  the  contour 
interval  is  10  units. 
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Figure  10.  Specific  Volume  Anomaly  (T),  July,  10  im  6f,  Units  arc  10  Ag 


and  the  contour  interval  is  10  units. 
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I0-*  kg-\  and  the  contour  interval  is  1  unit  (dashed  contour  interval 


is  0.5  units  from  solid  contours). 


500  m,  3t  (Figure  13)  and  Ss  (Figure  14)  fields  arc  virtually  mirror  images  of  each  other. 
At  that  depth,  ranges  of  40  to  70  are  typical  of  Os,  as  compared  to  45  to  75  for  Sj-.  At 
an  intermediate  depth  of  100  m,  ranges  from  65  to  190  and  70  to  150  correspond  to  Sr 
(Figure  15)  and  dj  (Figure  16),  respectively,  while  ranges  of  60  to  135  and  40  to  85  are 
representative  of  Sr  (Figure  17)  and  Ss  (Figure  18)  at  200  m.  The  salinity  signature  of 
the  Columbia  River  is  practically  gone  at  100  m.  In  fact,  in  the  region  north  of  43°  N, 
from  100  m  down,  gradients  in  both  the  Sr  and  Ss  fields  are  with  one  exception  consid¬ 
erably  weaker  than  in  southern  areas.  The  exception  is  a  5^  maximum  at  200  m  during 
October  in  the  data  rich  area  near  44°  N,  127°  W  (Figure  19).  This  une.xplained,  but 
possibly  real  (the  central  value  is  based  on  22  observations  although  surrounding  mean 
values  are  from  as  few  as  10)  mean  feature  supports  a  gradient  of  8.2  x  10-*  m’  kg-'  in 
one  degree  of  longitude,  which  is  comparable  to  the  strongest  gradients  an>■^vhere  in  the 
CCS  at  that  depth  and  much  stronger  than  any  other  in  the  northern  sector. 

At  100  m,  and  to  a  lesser  e.xtent  200  m,  a  meridional  tongue  of  5  r  with  its  axis 
just  offshore  is  apparent.  That  is,  there  is  a  pronounced  offshore  cold  core  adjacent  to 
the  California  coast.  Cyclonic  flow  around  this  feature  is  consistent  with  nearshore 
subsurface  poleward  flow  typical  of  the  C(  ^ 

C.  VARIANCE  AND  CORRELATION  OF  THE  SPECIFIC  VOLUME  ANOMALY 
CONSTITUENTS 

1.  Method 

Often  the  absolute  value  of  a  quantity  is  of  less  importance  than  the  nature  of 
the  temporal  or  spatial  derivatives  of  that  quantity.  Statistics  generated  from  a  signif¬ 
icant  number  of  observations  can  provide  information  concerning  those  differentials. 
Deviation  from  local  mean  values  can  cause  changes  in  local  gradient  or  prevent  change 
in  local  gradient  when  neighboring  values  change.  In  a  flow  supported  by  tight  gradi¬ 
ents,  small  meanders  result  in  relatively  large  changes  in  local  values.  The  actual  vari- 
ance2  of  a  field  is  a  descriptor  of  likely  deviation;  therefore,  in  a  data  block,  large 
variance  may  indicate  large  gradient  or  large  variability  in  the  local  gradient.  The  precise 
effect  on  instantaneous  gradient,  of  course,  depends  on  behavior  of,  and  correlation 
with,  neighboring  data  blocks.  For  example,  if  all  nearby  mean  values  change  in  tandem 
with  each  other,  the  gradient  would  remain  constant  and  "no  change"  in  gradient  would 
result  from  variance.  Without  thorough  knowledge  of  the  nature  of  the  field  of  interest, 

2  Whether  or  not  the  sample  variances  calculated  here  are  valid  representations  of  true  vari¬ 
ance  in  the  various  fields  is  a  relevant  question.  Section  IVFl  addresses  this  point  in  detail. 
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Figure  11  Specific  Volume  Anomaly  (T),  July,  500  m:  <5^.  Units  are 
10'*  m*  and  the  contour  interval  is  10  units  (dashed  contour  inter¬ 
val  is  5  units  from  solid  contours). 
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Figure  14.  Specific  Volume  Anomaly  (S),  July,  500  m:  Units  are 


10‘*  kg-*,  and  the  contour  interval  is  10  units  (dashed  contour  inter¬ 


val  is  5  units  from  solid  contours). 
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vSpecilk  Volume  Anomaly  (T),  July,  100  m:  <3,. 
10*'///’ A)?-',  anil  the  contour  interval  is  10  units. 


Units  arc 


K)-’  and  the  contour  interval  is  10  units. 
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Figure  17.  Specific  Volume  Aiioinnly  (T),  .luly,  200  m:  Or-  Units  arc 
10'* /«' /rg"',  and  the  contour  interval  is  IO  units. 
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10-*  m^kg-',  and  the  contour  interval  is  10  units. 
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Figure  19.  Specific  Volume  Anomaly  (S),  Oclolier,  2011  m:  6^  Units  arc 


and  titc  contour  ijitcrval  is  10  units. 
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it  is  impossible  to  quantify  gradient  from  sample  variance  calculation,  and  the  fact  that 
two  data  blocks  have  identical  variance  does  not  mean  that  they  have  identical  gradient. 
However,  determination  of  the  contribution  of  various  constituents  to  the  total  variance 
provides  strong  indication  of  the  contribution  of  those  constituents  to  the  development 
of  the  field  of  interest.  It  is  with  this  caveat  that  the  following  analysis  technique  is  used. 

Variance  of  the  specific  volume  anomaly  field  results  from  the  density  changes 
that  can  be  key  to  momentum  balances.  To  evaluate  the  importance  of  saline  variance 
to  specific  volume  anomaly  variance,  (variance  in  5s)  should  be  compared  to  a] 
(variance  in  5).  Because  equation  (3)  is  linear,  it  is  tempting  to  try  to  represent  the 
variance  of  5  as  the  sum  of  the  variances  of  the  constituents.  If  this  were  possible, 
comparison  of  ffjj.  (variance  in  dr)  to  aj,  for  example,  would  show  at  a  glance  how  much 
of  the  total  variance  is  due  to  temperature  variance  alone.  That  is,  whenever  the  ratio 
(7y<7j  remains  constant  and  large,  approaching  one  as  an  ideal  limit,  salinity  variance 
could  be  deemed  unimportant.  Conversely  if  varies  significantly  over  time  or 
space  intervals  of  interest,  salinity  variance  could  be  considered  important.  However, 
because  the  constituents  of  specific  volume  anomaly  are  not  independent  of  each  other, 
this  approach  is  not  valid.  The  covariances  play  an  essential  role  in  determining  the 
variance  of  the  specific  volume  anomaly  field.  Therefore,  to  evaluate  a  term  on  the  right 
hand  side  of  equation  (3)  as  a  potential  estimator  of  <5,  in  addition  to  variance  compar¬ 
ison,  it  is  necessary  to  determine  how  that  term  correlates  with  5.  If,  for  example,  dr  is 
highly  correlated  with  5,  then  most  of  the  variance  in  5  can  be  accounted  for  by  a  linear 
relationship  with  Sp  If  also  the  variance  of  <5r  is  of  similar  magnitude  to  that  of  5,  then 
5r  can  be  considered  a  good  estimator.  After  all,  if  5  changes  when  5r  changes,  in  similar 
amounts  and  directions,  then  5^  would  be  expected  to  model  d  wrell.  If  that  is  the  case, 
a  similar  look  at  the  behavior  of  5s  will  probably  show  it  contributes  little  to  total  vari¬ 
ance,  and  salinity  might  be  considered  unimportant  to  the  development  of  the  specific 
volume  field. 

The  correlation  coefficient  is  useful  as  an  indicator  of  the  degree  of  correlation 
but  more  useful  as  an  indicator  of  the  sign  of  the  correlation.  The  coefficient  of  deter¬ 
mination,  on  the  other  hand,  tells  exactly  how  much  of  the  variation  in  one  variable  is 
a  result  of  a  linear  relationship  with  the  other  (e.g.,  Walpole  and  Myers,  1985). 

These  ideas  can  be  used  to  evaluate  the  validity  of  the  use  of  temperature  alone 
to  determine  density.  Wherever  is  highly  correlated  with  5,  and  has  variance  of  sim¬ 
ilar  magnitude,  modeling  <5  with  dr  (or  p  with  p(T))  should  give  acceptable  results.  If 
5r  is  highly  correlated  and  its  contribution  to  total  variance  nearly  constant,  even  if  not 
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precisely  the  same  magnitude  as  total  variance,  a.  Sr  density  model  might  still  be  ac¬ 
ceptable  because  it  will  predict  a  representative  gradient  that  will  dificr  from  the  true 
gradient  by  a  nearly  constant  factor. 

2.  Results 

Variances  of  Sr,  and  Ss  were  computed  for  each  level  of  each  data  block,  as 
were  the  covariances,  (covariance  between  Sr  and  S)  and  (covariance  between 
Ss  and  S  ).  Correlation  coefilcients  between  S  and  Sr  and  S  and  Ss  were  also 
computed: 


V 

where  x  represents  either  Ss  or  Sp  In  addition,  coefficients  of  determination,  and 
(as  used  in  Figures  22-29)  were  calculated. 

Maps  of  determination  coefficients  (given  the*  sign  of  the  correlation  coefficient) 
for  all  months  and  the  depths  10,  100,  200  and  500  m  were  evaluated  (see  note  1). 
Multiplication  of  the  variance  in  5r  ,  by  r|^,  determines  how  much  of  the  total  variance 
in  S  is  due  to  variance  in  Sp  Below  about  100  m,  variance  in  the  S  field  is  very  small, 
too  small  to  allow  confidence  of  statistical  significance.  Analysis  of  Sr  and  5^  corre¬ 
lations  with  S  is  therefore  limited  to  the  upper  100  m.  Generally  Sr  (Figure  20, 
Figure  21)  is  more  highly  correlated  with  S  than  is  Ss-  Correlation  between  Ss  and  S 
(Figure  22,  Figure  23)  is  usually  quite  weak  on  the  surface  and  stronger,  although  not 
as  large  as  ,  at  100  m.  These  statements  are  true  only  in  a  very  general  sense.  For 
example,  frequently  in  both  fields  very  high  or  low  correlations  exist  in  the  immediate 
vicinity  of  each  other.  Another  exception  is  that,  as  expected,  immediately  offshore  of 
a  freshwater  source,  Sr  correlations  are  occasionally  weak,  or  nonexistent.  Usually  when 
the  Sr  correlation  is  weak,  the  Ss  correlation  is  very  high.  For  example  at  10  m  in  Feb¬ 
ruary  at  the  mouth  of  San  Francisco  Bay,  rj,,.  (Figure  24)  is  0.0  while  (Figure  25)  is 
1.0.  Both  correlations  can  be  very  high,  predictably,  for  instance,  in  the  Columbia  River 
plume,  where  warm,  fresh  water  results  in  high  specific  volume  anomaly.  On  the  other 
hand,  the  opposite  situation  exists  well  offshore  of  Baja  California,  near  Guadalupe  Is¬ 
land.  There,  in  late  winter  and  spring,  very  low  Sr  correlations  (Figure  26)  coexist  with 
very  low  values  (Figure  27).  Variance  in  both  Ss  and  dr  is  very  high  in  relation  to 
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5,  but  the  two  are  so  strongly  in  opposition  that  neither  correlates  well  with  the  total 
variance. 

Throughout  the  year,  within  a  two  or  three  degree  radius  of  31°  X,  119°  W  in 
the  upper  100  m,  r^j.  (Figure  28)  is  consistently  near  0.9  and  the  variance  in  Sj- 
(Figure  29)  accounts  for  a  reasonably  consistent  part  of  the  variance  in  d.  Perhaps  in 
this  area,  the  assumption  of  constant  salinity  would  introduce  minimal  error  in  a  density 
model. 


D.  MEAN  DIFFERENCES  OF  THE  SPECIFIC  VOLUME  ANOMALY 
CONSTITUENTS 
1.  Method 

If  an  observed  quantity,  .r,  is  modeled  by  a  predictor  x',  then  {x  —  x')  can  be 
considered  the  error  associated  with  a  single  observation.  If  that  error  is  squared  to 
eliminate  sign  bias  and  the  sum  of  such  squared  error  averaged,  the  square  root  of  the 
result  can  be  used  as  an  indicator  of  the  expected  deviation  and  is  termed  the  Root  Mean 
Square  Error  (RMSE).  The  RMSE  normalized  by  the  mean  value  of  then  provides 
an  estimate  of  the  magnitude  of  the  expected  error  in  terms  of  percent  of  the  "correct" 
value.  If  d  is  the  observed  quantity  and  the  predictor: 


—  RMSE„0,.,„^/. 


(13) 


If  for  example  x  is  T,  the  result,  RMSE(T),  provides  an  indication  of  how  well  Sr  ap¬ 
proximates  d  in  a  particular  data  block. 

2.  Results 

X’ormalized  RMSE  values  were  calculated  for  Sj-  and  Ss-  As  discussed  in  section 
IVB,  both  Sj-  and  ^5  are  everywhere  greater  than  zero  and  each  greater  than  the  absolute 
value  of  the  sum  of  the  remaining  specific  volume  anomaly  constituents.  Therefore  Sr 
and  Ss  are  necessarily  less  than  S,  and  the  sign  of  the  expected  error  is  known  (the  "root" 
and  "squared"  point  of  the  RMSE  technique  is  moot,  and  R.MSE  is  equivalent  to  mean 
dilTerence).  The  specific  volume  anomaly  will  always  be  underestimated. 

As  discussed  earlier,  it  is  the  gradient  in  the  error  field  that  is  likely  more  sig¬ 
nificant  than  magnitude  of  a  single  error,  all,  use  of  the  actual  mean  salinity  rather 
than  standard  35  psu  can  remove  a  constant  offset.  The  plots  of  R.MSE(T)  provide  a 
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Figure  20.  Specific  Volume  Cocllicicnt  of  Dcfenulnaliou  (  I ),  July,  10  m: 
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Figure  26.  Specific  VoUinie  Cocfricieiit  of  Dciciiiiin.’itiou  (T),  Maicit,  10  m: 
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Figure  27.  Specific  Volume  Coefficient  of  Determination  (S),  March,  iO  m: 
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Figure  28.  Specific  Volume  Coefficient  of  Determination  (T),  January,  10  m; 


Figure  29.  Variance  in  Steric  Anomaly  from  Temperature,  January,  10  m: 
>'Ur  100. 


clear  large  scale  picture  of  the  price  of  neglecting  salinity  variation  in  specific  volume 
anomaly  calculations  (e.g.,  Figure  30).  Near  the  surface  that  neglect  would  result  in  a 
40  to  50  percent  underestimation  in  the  north-south  gradient,  and  a  10  percent  under¬ 
estimate  west  to  east,  with  the  largest  errors  tending  to  occur  in  winter  months.  The 
meridional  gradient  persists  but  decreases  to  a  year-round  20  percent  at  500  m 
(Figure  31). 

Because  «5j-  and  ds  are  by  far  the  greatest  contributors  to  total  o,  data  block 
values  of  RMSE(S)  added  to  the  corresponding  values  for  dr  usually  are  nearly  equal  to 
one  (compare  Figure  32  to  Figure  30).  Nevertheless,  ds  mean  difference  plots  are  il¬ 
lustrative  of  the  part  ds  plays  in  determining  total  steric  anomaly.  The  large  scale  op¬ 
position  of  gradients  in  and  dr,  in  virtually  every  direction,  is  clear  from  analysis  of 
these  distributions.  At  500  m,  in  the  southern  section  of  the  CCS  there  is  a  pronounced 
increase  in  the  zonal  gradient  in  the  5^  RMSE  field  (Figure  33),  salinity  making  a  larger 
contribution  to  d  offshore  than  nearshore.  Some  smaller  scale  complexities  are  also 
discernible.  For  instance,  near  the  surface,  across  the  mouth  of  San  Francisco  Bay  in 
September  (Figure  34),  there  is  a  twenty  percent  change  in  the  dj  contribution  to  d,  ac¬ 
companied  by  a  very  flat  gradient  in  RMSE  for  dr-  This  phenomenon  is  illustrative  of 
the  local  nature  of  saline  variability,  and  temperature-salinity  interaction  in  the  CCS. 

An  interesting  result  of  this  analysis  is  that  there  are  large  areas,  virtually  any 
six  degree  latitude  by  four  degree  longitude  block,  where  the  assumption  of  no  salinity 
variation  would  apparently  result  in  less  than  a  ten  percent  error  in  the  resultant  specific 
volume  gradient.  That  is  particularly  true  if  areas  immediately  adjacent  to  the  coast  are 
excluded. 

E.  DYNAMIC  HEIGHT 
1.  Method 

Dynamic  height,  O,  is  a  measure  of  how  the  density  of  a  column  of  water  differs 
from  the  density  of  a  standard  column  of  the  same  depth  and  is  calculated  by  integrating 
the  specific  volume  anomaly  over  depth  from  an  assumed  layer  of  no  motion  to  a  depth 
of  interest.  It  is  the  horizontal  gradient  in  dynamic  height  that  determines  the 
geostrophic  velocity  field.  If  the  contributions  to  that  gradient  from  salinity  are  signif¬ 
icant,  they  should  be  considered  when  momentum  calculations  arc  performed.  If  the 
mean  differences  (5  —  Or),  where  (  )  denotes  the  mean,  in  columns  are  integrated  with 
respect  to  pressure,  the  resulting  horizontal  distribution  is  a  good  measure  of  the  effect 
of  salinity  on  the  geostrophic  velocity  field.  That  is,  contours  of  -  Oy,,  =  ^  define 
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Figure  31.  Specific  Volume  Anomaly  RMSE  (T),  July,  500  m. 
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Figure  32.  Specific  Volume  Anomaly  RMSE  (S),  .fuly,  10  m. 
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Figure  33.  Speciiic  Voltiinc  Anomaly  RMSE  (S),  .Inly,  500  m. 
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mean  geostrophic  flow  that  will  not  be  accounted  for  if  5  r  is  used  as  an  estimator  of  5, 
Therefore,  if  gradients  in  the  difference  field  are  large,  salinity  is  important. 

2.  Results 

Plots  of  monthly  mean  values  of  O,  Oj.,,,  and  O',  obtained  by  integrating  from 
an  assumed  layer  of  no  motion  at  500  m,  to  200  m  (200/500),  100  m  (100/500),  and  the 
surface  (0/500)  were  plotted  and  evaluated.  The  O  results,  where  they  overlap,  are  in 
ve:y  close  agreement  with  those  of  Wyllie  (1966).  Standard  deviations  are  in  virtually 
every  case  considerably  more  than  an  order  of  magnitude  less  than  mean  values.  Even 
so,  this  small  variance  still  indicates  that  contour  locations  are  doubtful  by  more  than 
the  width  of  a  contour  inter\'al  on  these  maps.  Therefore,  any  conclusions  drawn  con¬ 
cerning  small  scale  features  are  suspect.  It  should  also  be  remembered  that  these  are 
maps  of  mean  fields,  representing  climatology  rather  than  instantaneous  flows. 

Examination  of  the  0/500  dynamic  height  field  of  the  CCS  (e.g..  Figure  35  and 
Figure  36),  leads  to  division  of  the  area  into  three  regions.  North  of  43®  N  the  diver¬ 
gence  of  the  North  Pacific  Current  and  freshwater  discharges  from  the  Columbia  River 
and  Strait  of  Juan  de  Fuca  result  in  disorganized  mean  geostrophic  gradients.  A  split 
of  the  flow  toward  the  north  and  south  is  superimposed  on  weak  circulations  around 
warm,  freshwater  plumes.  The  area  along  the  U.S.  coast  from  30°  to  38°  N  is  dominated 
by  cyclonic  geostrophic  circulation  inshore,  and  flow  generally  to  the  southeast  offshore. 
The  strength  of  the  southward  flow  and  degree  of  cyclonic  turning  onshore  is  seasonal. 
Southward  flow  is  strongest  in  summer  and  weakest  in  winter.  Conversely,  onshore 
cyclonic  turning  is  least  in  summer,  greatest  in  winter.  There  is  another  cyclonic  center 
just  offshore  of  Baja  California  near  25°  N.  Again  the  cyclonic  departure  from  a  gen¬ 
erally  southeastward  flow  is  seasonal  as  described  above.  For  each  month  the  dynamic 
height  difference  fields  (O')  (Figure  37)  show  inclusion  of  horizontal  salinity  gradients 
to  be  responsible  for  a  nearshore,  equatorward  geostrophic  component  from  San 
Francisco  to  San  Diego  and  a  southwestward  flow  along  the  coast  of  Baja  California. 
This  predictable  geostrophic  pattern  is  due  to  the  persistent  salinity  pattern:  fresh  to 
saline,  offshore  toward  shore,  with  a  more  meridional,  fresh  to  salty  north  to  south, 
gradient  in  the  southern  section.  What  is  even  more  interesting  is  the  fact  that 
throughout  the  year,  the  magnitude  of  the  anomalous  flow  implied  by  the  <!>'  field  is  on 
the  same  order  as  flow  associated  with  VO  itself,  and  perhaps  greater  in  some  areas  than 
the  magnitude  of  the  flow  predicted  by  VOjj,^  (compare  figures  36-38). 

In  the  northern  region,  in  each  month,  VOj<,^  fails  to  account  for  a  persistent 
south-southwestward  component  from  the  mouth  of  the  Columbia  River  that  varies  in 
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Figure  35.  Dynamic  lleiglit,  0/500  m,  .laiuiary:  <I),  the  contour  interval  is  0.2  and 
units  arc  dynamic  decimeters. 
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Figure  37.  Dyn.iinic  Meiglit  Difference,  0/500  in,  .Inly:  <!>'  =  <I>  —  the  con¬ 
tour  interval  is  0.2  and  units  are  dynamic  decimeters. 


Figure  38.  Dynamic  lleiglit  (Salinity  35),  0/500  in,  .luiy:  tlic  contour  inter¬ 

val  is  0.2  and  units  arc  dynamic  decimeter.':. 
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intensity  with  the  magnitude  of  discharge  from  that  river.  However,  seaward  of  that 
southward  component  the  difference  field  is  very  slack.  That  is,  the  dynamic  height  field 
estimated  from  temperature  gradient  alone  verj'  nearly  accounts  for  the  geostrophic  flow 
pattern  in  the  region  where  surface  salinity  has  the  largest  annual  variance.  This  ap¬ 
parent  contradiction  is  explained  by  the  fact  that  the  horizontal  salinity  gradient  here 
decreases  rapidly  with  depth.  At  depth,  throughout  the  year,  there  is  virtually  no  sig¬ 
nificant  salinity  gradient  in  this  northern  area.  At  200  m  for  example,  over  a  distance 
where  mean  temperature  may  vary  by  0.4°  C,  salinity  changes  less  than  0.1  psu.  In  the 
southern  area,  by  contrast,  a  0.5°  C  temperature  change  is  accompanied  typically  by  a 
change  of  0.3  psu.  The  result  is  that  in  the  north,  the  depth  integrated  effect  of  salinity 
on  geostrophic  flow  only  approaches  that  of  temperature  in  the  immediate  vicinity  of  the 
sharpest  upper  level  salinity  gradient. 

Analyses  of  dynamic  height  and  geostrophic  flow  at  other  levels  {e,g„  200/500 
and  100/500  m,  Figures  39-44)  show  relative  gradients  very  similar  to  those  previously 
discussed.  VO'  is  of  similar  order  as  VO  and  VOjj,,  at  all  depths,  and  the  general  pattern 
of  the  diflerence  field  persists  with  depth.  North  of  northern  Baja  California,  the  differ¬ 
ence  field  shows  a  strong  offshore-onshore  gradient  indicative  of  higher  salinities  inshore 
in  the  poleward  flow.  This  result  is  consistent  with  the  idea  that  the  poleward  flow  can 
be  traced  by  its  high  salinity  {e.g.,  Lynn  and  Simpson,  1987).  That  is,  the  inshore 
poleward  undercurrent  may  carry  the  saline  signature  of  Equatorial  water  which  is  to  the 
south  of  the  CCS.  That  salinity  increases  significantly  toward  the  south  is  evident  from 
the  fact  that  south  of  about  34°  N  there  is  a  strong  alongshore  north  to  south  gradient 
in  the  difference  field.  The  southward  salinity  increase  is  also  apparent  in  the  salinity 
fields  (Figure  4).  The  100/500  and  200/500  m  dynamic  height  calculations  show  in¬ 
creased  cyclonic  turning  associated  with  the  large  scale  California  coastal  features  de¬ 
scribed  earlier.  That  is,  the  dense  cores  associated  with  what  look  like  "c"  shaped  troughs 
centered  at  the  California  coast  and  Baja  California  peninsula  on  the  0/500  m  plots  move 
offshore  with  depth  so  that  at  200  m,  and  possibly  even  at  100  m,  a  closed  geostrophic 
circulation  with  coastal  poleward  flow  exists.  Because  of  the  neglected  zonal  salinity 
gradient,  greatly  overestimates  the  strength  and  horizontal  extent  of  the  poleward 
flow. 

Due  to  the  absence  of  the  freshwater  runoff  below  the  surface,  the  divergence 
of  the  Nortli  Pacific  Current  is  relatively  clear  on  the  100/500  and  200;  500  m  plots.  This 
divergence  area  seems  to  be  farther  south  in  the  winter  than  summer,  as  noted  by 
Pickard  and  Emery  (1982).  Although  the  divergence  occurs  in  a  region  of  generally  poor 
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data  coverage,  the  seasonal  migration  seems  to  be  from  near  40®  N  in  the  winter  to 
45®  N  in  summer,  rather  than  the  45®  N  to  50®  N  range  reported  by  Pickard  and  Emery 
(1982). 

The  focus  of  this  thesis  is  large  scale;  however,  using  the  rule  of  thumb  that 
features  whose  wavelengths  exceed  2  grid  points  in  extent  may  be  resolved,  it  may  be 
possible  to  use  these  calculations  to  gain  insight  into  some  relatively  small  scale  features. 
In  Appendix  B,  discussion  of  some  smaller  scale  effects  of  salinity  on  dynamic  height  in 
the  CCS  is  provided. 

F.  PROPAGATION  OF  T-S  VARIANCE  THROUGH  THE  EOS 
1.  Method 

An  alternate  approach  to  examine  the  relative  importance  of  constituent  vari¬ 
ance  to  density  variance  is  to  propagate  the  variance  associated  with  each  of  the  inde¬ 
pendent  variables  through  the  complete  nonlinear  equation  of  state.  By  the  chain  rule, 
if  density  is  strictly  a  function  of  salinity,  temperature,  and  pressure  then: 


A  A  ,  A  ,  K 

dip  =  -7—  A/  +  -:r"  As  -I-  -T—  A/7, 
ct  cs  op 


(14) 


where  (A)  indicates  departure  from  mean  values  and  p  is  mean  density.  Because  all  ob¬ 
servations  arc  on  standard  depth  levels  and  pressure  changes  on  those  levels  are  small, 
even  in  comparison  to  error  associated  with  depth  determination,  Ap  can  be  neglected. 
Squaring  both  sides  of  (14),  summing  over  a  number  of  observations,  and  dividing  by 
the  number  of  observations  less  one  (e.g.,  Greenwalt  and  Shultz,  1962)  leads  to  : 


,  cp  v2  a  .  ,  ^2  2  ,  w  cp 


(15) 


where  represents  variance  (either  sample  or  population),  a  denotes  covariance,  and 
subscripts  represent  variables  to  which  the  statistic  applies.  As  demonstrated  by 
Lillibridge  (1989),  the  equation  of  state,  EOS(80),  a  polynomial,  is  relatively  easy  to 
differentiate.  Therefore  the  contribution  to  n\  from  temperature,  salinity,  and  their 
interaction  is  easily  determined  over  any  interval  of  interest.  This  technique  can  be  as 
readily  applied  to  specific  volume  or  any  of  the  specific  volume  anomaly  constituents 
and  is  desirable  because  it  employs  the  full  nonlinear  equation  of  state.  Insofar  as  the 
statistics  used  are  accurate  estimators,  the  equation  of  state  itself  accurate,  and  the 
pressure  levels  constant,  the  c.xprcssion  is  valid.  The  method  is  particularly  useful  when 
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Figure  41.  Dynamic  Height  Difference,  l()0/5()0  in,  .luly:  d)' =  <I)  —  the 
contour  interval  is  0.2  and  units  arc  dynamic  decimeters. 
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Figure  42.  Dynamic  Height,  200/500  m,  .Inly:  (!>,  tl)c  coiuour  interval  is  0.2  and 
units  arc  dynamic  decimeters. 
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Figure  43.  Dynamic  Height  (Salinity  35),  200/501)  ni,  July:  (D,,,,,  the  contour  in- 
terval  i$  0.2  and  units  are  dynamic  decimeters. 
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temperature  and  salinity  ranges  over  a  time*space  inter\’al  of  interest  are  reasonably  well 
known  and  a  quick  assessment  of  the  relative  importance  of  constituent  variances  de¬ 
sired.  It  also  is  computationally  quite  eflicient  because  the  equation  of  state  need  only 
be  invoked  once  to  determine  total  and  constituent  variance  over  a  sample  of  interest. 

Let  the  terms  on  the  right  hand  side  of  equation  (15)  be  represented  as: 


2  f  (’P  y,2  2 


op 

ds 


To: 


(16) 
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As  mentioned  earlier,  to  increase  the  number  of  observations  and  hence  statistical  va¬ 
lidity,  these  data  block  variance  estimates  were  computed  over  three-month  seasons, 
with  winter  defined  as  January,  Februaiys  .March. 

Potential  energy  {PE)  of  a  column  of  water  can  be  determined  by  integrating 
PE  =  pgz,  where  g  is  the  acceleration  due  to  gravity,  considered  constant,  and  z  is  height 
above  some  reference  level,  over  depth  (e.g.,  FofonolT,  1962).  It  follows  that  on  a  con¬ 
stant  depth  level  (or  very  thin  layer)  p  is  the  only  variable  in  PE,  and  the  horizontal 
density  field  can  be  used  to  represent  potential  energy.  Change  in  the  density  field,  either 
over  time  or  distance,  then  represents  change  in  the  potential  energy  field.  Equivalently, 
variance  from  the  mean  density  field  can  be  considered  as  an  indication  of  energetic  ac¬ 
tivity  or  potential.  In  general,  relatively  high  density  variance  should  correspond  to 
relatively  high  energetic  activity.  Therefore,  if  sample  variance  computed  over  some  time 
or  space  interval  is  an  accurate  estimator  of  the  true  density  variance,  sample  variance 
may  provide  indication  of  energetic  activities.  Dantzler  (1977)  used  similar  statistical 
techniques  applied  to  vertical  excursions  of  the  thermocline  to  calculate  North  Atlantic 
potential  energy  levels. 

2.  Results 

.Maps  of  oj  and  percent  contributions  of  salinity,  (aljol  x  100),  temperature, 
{aljcl  X  100),  and  covariance,  x  100),  were  evaluated.  Near  the  surface,  large 

density  variance  occurs  in  areas  of  maximum  fluctuation  in  precipitation,  heating,  or 
freshwater  influx  for  example.  Accordingly,  Figure  45  indicates  high  variance  near  the 
Columbia  River  outflow,  and  in  the  upwelling  region  along  the  coast  of  California.  As 
expected,  salinity  variance  is  dominant  in  the  Columbia  River  plume  (Figure  46).  Also 
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as  anticipated,  in  the  upwelling  region,  temperature  and  salinity  effects  (compare 
Figure  47  and  Figure  46)  are  positively  correlated  and  of  similar  relative  importance  to 
total  density  variance.  These  near-surface  observations  are  consistent  with  results  al¬ 
ready  presented  and  so  provide  evidence  that  this  technique  may  be  valid. 

Away  from  the  surface,  areas  of  highest  indicated  variance  may  be  dynamically 
active  regions  in  the  CCS.  Regions  of  maximum  might  correspond  to  the  mean  po¬ 
sitions  of  currents  or  eddies.  If  so,  equation  (15)  can  be  used  to  determine  the  relative 
contributions  of  temperature,  salinity,  and  their  interaction,  to  the  CCS  itself.  In  sup¬ 
port  of  this  idea,  the  locations  of  maximum  aj  seem  to  more  closely  correspond  to  CC 
and  CU  locations  as  depth  increases.  For  example  areas  with  trj  at  least  twice  as  high 
as  neighboring  data  blocks  at  200  m  (Figure  48)  correspond  reasonably  well  with  the 
areas  of  tight  gradient  in  O  (200/500  m)  apparent  in  Figure  42.  In  the  regions  of  high 
variance,  salinity,  with  variance  ratios  (oy/cj  typically  from  50  to  200,  seems  to  con¬ 
tribute  more  to  total  variance  than  temperature,  with  variance  ratios  {olJaD  in  the  high 
total  variance  region  of  5  to  150.  At  500  m,  salinity  contributes  virtually  all  of  the  density 
variance  (Figure  49)  in  a  region  that  resembles  the  characteristic  cyclonic  circulation  of 
the  CCS.  The  high  correlation  of  saline  contribution  to  density  variance  indicates  that 
the  density  gradient  is  largely  due  to  salinity  gradients  at  this  depth.  Figure  50  shows 
some  areas  at  500  m  with  total  density  variance  values  on  the  same  order  as  at  200  m 
and  these  areas  have  values  at  least  twice  as  high  as  neighboring  data  blocks.  The  lo¬ 
cation  of  these  regions  match  the  location  of  some  of  the  energetic  portions  of  the  CCS 
discussed  above.  These  observations  provide  evidence  that  the  density  gradient  at 
500  m  may  be  significant  and  is  largely  induced  by  salinity  gradients.  This  speculation 
is  consistent  with  that  of  Lynn  and  Simpson  (1989). 

If  quarterly  mean  values  of  salinity  and  temperature  are  substituted  for  t  and  s 
in  equation  (15),  and  statistics  are  generated  from  these  mean  values  by  summing  across 
the  length  or  width  of  the  CCS,  equation  (16)  will  represent  contributions  of  salinity, 
temperature,  and  their  covariance  to  the  total  density  variance  and  may  indicate  their 
importance  to  the  mean  density  gradient.  Such  calculations  were  performed  to  obtain 
statistics  cross-shore  at  each  latitude  that  had  at  least  five  seasonal  means  of  at  least  10 
observations,  and  alongshore,  across  data  blocks  equidistant  from  the  coast. 
Alongshore  lines  are  numbered  1  to  8  from  inshore  out,  and  cross-shore  lines  are  iden¬ 
tified  by  latitude  (Figure  51).  The  results  at  the  surface  and  100  m  v.'cre  plotted  as 
histograms.  Inshore,  alongshore  lines  I  (Figure  52)  and  2  (Figure  53)  show  large  year- 
round  total  variance  at  the  surface.  In  the  northern  sector  where  line  1  crosses  the 
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Figure  47.  Tciuperadirc  Coni rilnii ion  to  Tot.il  Density  Variance,  It)  in, 
Suinincr:  x  100. 
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Figure  50.  Tot.il  density  Variance,  500  in,  Suiniiicr:  Units  arc  {kg  iir^y  x  |0^ 
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Columbia  River  outflow  region,  the  greatest  surface  contribution  is  clearly  from  salinity. 
On  this  line  even  at  100  m,  salinity  variance  dominates  alongshore  except  in  winter  when 
temperature  variance  is  apparently  slightly  more  important.  Line  2,  which  is  influenced 
by  runoff  the  length  of  the  California  coast,  as  well  as  by  Columbia  River  outflow,  shows 
surface  salinity  as  most  important  in  fall,  and  a  strong  saline  signature  year  round.  At 
100  m,  temperature  is  clearly  dominant  in  the  low  runoff,  winter  season,  and  at  least  of 
equal  importance  the  rest  of  the  year.  Northern  cross-shore  sections  (only  sections  at 
46°,  47°,  49°  N  in  the  autumn,  and  26°-35°  N  met  the  minimum  number  of  observations 
criterion)  show  salinity  variance  as  key  during  the  fall  (e.g..  Figure  54). 

Alongshore,  the  covariance  is  always  negative,  as  must  be  expected  when  the 
large  scale  opposing  meridional  trends  in  temperature  and  salinity  are  considered. 
Cross-shore,  except  in  the  northern  area  where  warm,  fresh  Columbia  River  outflow  re¬ 
sults  in  strong  positive  interactions,  covariance  tends  to  be  less  important  and  of  either 
sign,  reflecting  the  variability  of  the  nearshore  waters  (e.g..  Figure  55). 

Cross-shore  variance  at  100  m  is  vciy'  small  in  the  north,  generally  an  order  of 
magnitude  less  than  the  O.l  to  0.5  {kgm-^y  typical  at  the  surface.  In  the  south,  cross¬ 
shore  variance  at  100  m  is  of  comparable  magnitude  to  that  at  the  surface. 

Histograms  generated  for  areas  where  relative  contributions  to  density  variance 
of  temperature,  salinity,  and  their  interaction  are  the  same  will  have  similarly  shaped 
histograms.  Therefore,  it  is  reasonable  to  expect  that,  in  the  same  water  mass,  areas 
undergoing  similar  energetic  processes  will  have  similarly  shaped  histograms.  This 
reasoning  implies  that  if  the  shape  of  the  histogram  for  a  particular  section  is  the  same 
at  the  surface  and  100  m,  the  two  levels  may  well  be  undergoing  the  same  dynamics 
and/or  thermodynamics;  that  is,  they  are  dynamically  coupled.  In  the  autumn,  the 
graphs  for  areas  south  of  31°  X  at  both  depths  arc  vcr>'  similar  in  shape  (e.g.. 
Figure  56),  indicating  such  coupling  may  exist.  Perhaps  surface  cooling  which  occurs 
during  the  months  of  October  through  December  results  in  deeper  mixed  layers  south 
of  31°  X. 

In  this  southern  region,  except  when  total  variance  is  very-  low,  temperature 
provides  the  largest  contribution  to  cross-shore  gradients.  There  is  a  tendency  for 
greatest  temperature  variance,  and  therefore  total  variance,  to  occur  in  fall,  further  evi¬ 
dence  that  surface  cooling  is  greatest  during  that  period. 

Alongshore  lines  also  show  a  noticeable,  if  less  pronounced,  increase  in  variance 
during  the  autumn.  During  that  season,  the  alongshore  sections  seaward  of  line  2  show 
the  surface  and  100  m  layers  to  be  well  coupled  the  length  of  the  CCS  (e.g..  Figure  57). 
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Figure  51.  Alongshore  and  Cross-Shore  Density  Variance  Contribution:  Con¬ 
tributions  to  total  density  variance  by  temperature,  salinity,  and  their 
interaction  were  calculated  along  the  lines  shown. 
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Figure  52.  Alongshore  Line  I,  Conirihufions  to  Density  Variance;  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
Ill.  Units  of  density  variance  here  arc  (Agur^y. 
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m.  Units  of  density  variance  here  arc  {kg  »r’)^ 


76 


Figure  54.  Cross-shore  Line  46,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  licrc  by  season  at  the  surface  and  100 
ni.  Units  of  density  variance  here  arc  {kgnr^y. 
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Figure  55.  Cross-shore  Line  26,  Confributions  io  Density  Variance:  Coniribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surrace  and  100 
in.  Units  of  density  variance  here  arc  {kg  nr^y. 
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Again  a  possible  interpretation  is  that  surface  cooling  in  the  fall  results  in  deeper  mixing. 
If  this  interpretation  is  valid  then  mixing  occurs  offshore  for  the  length  of  the  CCS. 

It  may  be  argued  that  over  any  areas  or  periods  where  the  relative  contributions 
to  total  variance  remain  constant,  a  density  model  based  on  temperature  alone  can  ac¬ 
curately  represent  gradients,  erring  only  by  a  constant  factor.  In  that  case,  intersections 
of  a  group  of  adjacent  alongshore  sections  sharing  a  characteristic  histogram  shape,  with 
a  group  of  cross-shore  sections  that  have  a  common  histogram  shape,  will  highlight  such 
regions.  Subject  to  individual  interpretation,  some  seasonal  combinations  of  lines  2-5 
and  26-31  may  meet  the  criteria.  For  example,  the  autumn  histograms  for  lines  3 
(Figure  58)  and  4  (Figure  59)  are  nearly  the  same  shape,  as  are  the  autumn  graphs  of 
lines  27,  28,  and  29  (Figure  60,  Figure  56  and  Figure  61,  respectively).  The  implication 
here  is  that  during  October  through  December,  at  least  in  the  upper  100  m,  within  an 
area  2‘’-4°  from  shore,  between  27°  and  30°X,  a  temperature-only  density  model  would 
overestimate  the  alongshore  gradient  by  a  constant  small  factor  and  underestimate  the 
cross-shore  gradient  by  a  similar  amount.  Again  the  indication  is  that  offshore,  in  select 
three  or  four  degree  blocks,  salinity  variance  may  not  be  essential  to  an  accurate  model 
of  the  dynamic  field.  (Using  this  technique  to  generate  shorter  alongshore  cuts,  of 
lengths  similar  to  the  cross-shore  cuts  for  example,  might  result  in  the  identification  of 
more  of  these  areas.) 


Figure  56.  Cross-shore  Line  2.S,  Coiilriluitioiis  (o  Dciisiiy  Vari.incc:  ('ontribiition 
to  lotol  dcnsily  variance  (labeled  T)  by  teinperaUiic  (1),  salinity  (s),  and 
their  interaction  (l/s)  are  depicted  here  by  season  at  the  suiiace  aiul  lOO 
in.  Units  ol'den.sity  variance  here  arc 
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Figure  57.  Alongshore  Line  6,  Contrihutions  (o  Density  Variance:  Contribution 
to  total  density  variance  (labeled  ’I  )  by  temperature  (t).  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  licre  by  season  at  tiic  surCacc  and  100 
in.  Units  of  density  variance  here  arc  {k}>  nr^y. 
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Figure  58.  Alongshore  Line  3,  Contributions  to  Density  Variance:  Contribution 


to  total  density  variance  (labeled  T)  by  temiicralurc  (l),  salinity  (s),  and 
their  interaction  (t/s)  are  depicted  here  by  season  at  the  surface  and  100 


ni.  Units  of  density  variance  here  are  (%  nr^y. 
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Figure  59.  Alongshore  Line  4,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  are  depicted  iicrc  by  sc3soi\  at  the  suiface  and  100 
m,  Units  of  density  variance  iicre  arc  {kg  nr^y. 


Figure  60.  Cross-shore  Line  27,  ContriluOions  to  Density  Varj.ince!  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
m.  Units  of  density  variance  here  arc  {kg  m-y. 
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Figure  61.  Cross-shore  Line  29,  Contrihutiotis  to  Density  Variance:  Contribution 


to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 


their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 


ni.  Units  of  density  variance  here  are  {kgur^. 
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V.  SUMMARY  AND  CONCLUSIONS 


Salinity  has  a  profound  impact  on  the  CCS,  with  structure  that  is  complex  and 
variable.  In  the  region  of  the  Columbia  River's  outflow,  surface  salinity  may  var>'  an 
entire  psu  in  100  km,  or  1  psu  over  the  course  of  a  year;  yet  in  this  region  the  depth  in¬ 
tegrated  importance  of  the  salinity  gradient  is  less  important  than  anj'where  else  in  the 
CCS.  Off  the  coast  of  Baja  California,  where  thermal  expansion  coefficients  are  high, 
saline  contraction  coefficients  low,  and  surface  salinity  gradients  have  minimal  impact 
on  the  density  field,  the  importance  of  salinity  variability  100  m  below  the  surface  is  of¬ 
ten  equal  to  the  importance  of  temperature  variance.  The  mean  salinity  field  has  been 
shown  to  have  a  pronounced  effect  on  the  large-scale  geostrophic  circulation  of  the  CCS. 
On  the  other  hand,  there*seem  to  be  regions  on  the  order  of  3°  squares  where  the  as¬ 
sumption  of  constant  salinity  might  not  cause  significant  errors  in  the  estimation  of 
density  gradients.  Relative  importance  of  temperature  and  salinity  often  changes  sea¬ 
sonally.  Areas  where  mean  gradients  in  temperature  and  salinity  have  positive  corre¬ 
lation  in  terms  of  their  effect  on  density  border  areas  of  negative  covariance  of  the  T-S 
gradients.  Despite  the  existence  of  correlation  in  mean  gradients,  there  seems  to  be  little 
correlation  between  simultaneously  observed  values  of  temperature  and  salinity. 

Based  on  these  observations,  it  is  reasonable  to  conclude  that  careful  work  con¬ 
cerning  the  CCS  should  include  haiine  description.  Descriptively  and  dynamically, 
salinity  is  essential  to  accurately  characterize  the  large-scale  structure  of  the  CCS.  Evi¬ 
dence  has  been  developed  in  this  study  that  indicates  salinity  may  play  significant  roles 
in  many  smaller  scale  phenomenon  as  well.  There  is  some  indication  that  there  may  be 
areas  where  salinity  variation  may  be  neglected  in  process-oriented  models,  but  regions 
to  be  so  modeled  must  be  carefully  selected. 

A.  STERIC  ANOMALY 

A  comparison  of  the  constituents  of  the  steric  anomaly  due  to  salinity  and  temper¬ 
ature  showed  that  the  gradients  in  dj  and  ds  oppose  generally,  major  exceptions  being 
on  the  surface  off  the  coast  of  Oregon  where  the  warm,  fresh  discharge  of  the  Columbia 
River  leads  to  large  values  of  dy  and  ds,  and  off  the  coast  of  northern  California  in 
summer,  when  upwelling  of  cold,  salty  water  leads  to  low  values  for  both  dr  and  5^. 
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In  areas  where  5^  is  highly  correlated  with  6,  and  has  similar  variance,  Sj-  might  be 
considered  an  adequate  estimator  of  3.  A  small,  shallow  region  with  these  characteristics 
was  identified  off  the  coast  of  northern  Baja  California. 

Using  the  idea  of  RMSE  to  calculate  mean  differences  from  the  steric  anomaly,  S, 
of  and  St,  the  contributions  to  mean  gradient  in  5  associated  with  Ss  and  Sr  were 
compared.  On  the  large  scale,  is  quite  important  and  certainly  should  not  be  ignored. 
On  the  other  hand,  the  contribution  of  salinity  is  fairly  evenly  distributed  and  so,  as  the 
amount  of  area  considered  decreases,  the  relative  magnitude  of  error  in  gradient  resulting 
from  the  assumption  of  constant  salinity  also  decreases.  For  example,  assumption  of  a 
constant  value  of  salinity  for  the  entire  CCS  in  July  would  result  in  an  approximate 
overestimate  of  40  percent  in  the  alongshore  south  to  north  mean  S  gradient.  However, 
if  only  the  area  adjacent  to  the  California  coast  is  considered,  a  10  percent  error  results. 

B.  DYNAMIC  HEIGHT 

Calculating  dynamic  height,  <I>,  by  integrating  5  over  depth,  and  comparing  it  to 
dJjj,,  ,  calculated  by  integrating  St  ,  provided  insight  concerning  the  importance  of 
salinity  variation  to  the  geostrophic  flow  field.  At  all  depths,  the  mean  variability  of 
salinity  in  the  CCS  is  responsible  for  a  significant  equatorward  component  along  the 
coast  of  California  and  a  strong  offshore  component  adjacent  to  Baja.  Additionally 
mean  local  salinity  giadicnts  can  play  a  major  role  in  the  structure  of  small  scale 
features. 

C.  VARIANCE  IN  THE  DENSITY  FIELD 

By  calculating  temperature  and  salinity  variances  and  covariances  and  differentiating 
the  equation  of  state,  contributions  to  local  variation  in  density  due  to  salinity,  temper¬ 
ature,  and  their  interaction  were  computed  directly.  Temperature  variation  provided  the 
major  contribution  year-round  on  the  surface.  Salinity's  importance  generally  increased 
with  depth  so  that  at  500  m,  salinity  variance  might  actually  be  most  significant.  Sub¬ 
stantial  variance  in  the  density  field  existed  even  at  500  m  in  locations  that  correspond 
to  the  CC  and  CU,  indicating  dynamic  activity  was  occurring  at  500  m,  and  perhaps 
that,  as  Lynn  and  Simpson  (1989)  have  suggested,  the  CCS  can  extend  deeper  than 
500  m.  and  that  at  deeper  depths  it  may  be  traceable  by  salinity  variance. 

By  applying  this  technique  to  the  mean  temperature  and  salinity  field,  the  relative 
contribution  to  longshore  and  cross-shore  gradient  was  determined.  Nearshore,  or 
elsewhere  in  the  vicinity  of  river  drainage  (Figure  52).  variance  is  generally  high  and 
saline  contribution  considerable.  There  is  a  tendency  for  high  variance  to  occur  in  t.ic 


87 


autumn  (defined  here  as  October  through  December,  Figure  57)  due  apparently  to  sur* 
face  cooling.  Offshore  this  cooling  and  associated  mixing  tend  to  decouple  the  surface 
layer  and  allow  greater  density  variation  at  deeper  depths,  at  least  100  m. 
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VI.  RECOMMENDATIONS 


A.  PROFILE 

The  information  developed  in  this  research  should  bc  displayed  in  vertical  sections. 
The  organizational  techniques  employed  here  would  allow  alongshore,  cross-shore,  or 
diagonal  lines  to  be  displayed  easily  in  vertical  cross  section.  Study  of  these  sections 
would  allow  fuller  appreciation  of  the  three  dimensional  mteraction  of  temperature  and 
salinity  in  the  CCS.  Perhaps  sections  of  the  contribution  of  salinity  to  density  variance 
could  be  used  wth  the  ideas  of  Lynn  and  Simpson  (1989)  concerning  "spiciness"  to  gain 
a  fuller  understanding  of  the  water  mass  structure  of  the  CCS. 

B.  FINER  RESOLUTION 

This  analysis  has  been  carried  out  at  a  1“  spatial  by  1  month  temporal  scale,  and  so 
allows  only  large-scale  resolution.  A  much  finer  scale  approach,  perhaps  r;2®  by  2 
weeks,  might  be  possible  in  certain  regions.  Using  similar  processing  techniques  and 
finer  resolution,  if  eveiy  depth  level  were  analyzed,  rather  than  the  few  depths  empha¬ 
sized  here,  the  T-S  interplay  associated  with  small  scale  phenomena  such  as  coastal 
upwelling  might  be  characterized. 

C.  MODELING  CONSIDERATIONS 

Speculation  has  been  offered  here  concerning  the  validity  of  the  assumption  of 
constant  salinity  for  the  purpose  of  dynamic  modeling.  Ideally,  the  mean  salinity  data 
processed  for  this  study  should  be  incorporated  into  a  CCS  model  run  and  results  com¬ 
pared  to  similar  runs  in  which  salinity  was  held  constant.  On  the  other  hand,  indications 
of  the  effects  of  salinity  might  be  achieved  with  less  effort  by  using  a  temperature-driven 
linearization  of  density  that  simulates  saline  variability.  For  example,  increasing  the 
magnitude  of  the  thermal  expansion  coefficient,  a,  in  an  upwelling  process  model  may 
capture  the  effect  of  a  positive  T-S  covariance  on  density.  Entering  different  values  of 
a  for  different  regions  and  depths  (higher  values  where  covariance  is  positive,  lower  val¬ 
ues  where  negative)  would  simulate  salinity's  variation  over  regions  where  salinity  or 
temperature  effects  change  sign.  Use  of  a  constant  value  of  a  over  a  range  of  temper¬ 
atures,  in  addition  to  ignoring  saline  contribution,  can  superimpose  a  false  smoothing 
tendency,  resulting  in  an  underestimation  of  the  magnitude  of  positive  Ap  in  cold  water, 
and  negative  Ap  in  warm  water.  Forcing  a  to  change  more  or  less  uniformly 
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meridionally  but  at  a  higher  rate  than  warranted  by  the  mean  temperature  gradient 
would  allow  study  of  the  combined  effects  of  meridional  saline  variation  and  smoothing. 

Because  salinity  measurements  arc  more  expensive  to  make  and  more  difficult  to 
obtain  than  temperature  observations,  it  may  be  desirable  to  estimate  salinity  from  a 
known  temperature-salinity  (T-S)  correlation  (Stommel,  1947).  Emeiy  (1975)  has  shov^m 
that  the  accuracy  of  dynamic  height  calculations  made  from  such  derived  salinities  can, 
at  least  in  some  cases,  be  comparable  to  the  accuracy  of  determinations  made  from  ob¬ 
served  salinities.  In  these  cases,  the  salinity  field  in  general  can  be  accurately  inferred 
from  expendable  bathythermograph  observations  alone,  as  Rienecker  et  al  (1985)  have 
attempted. 

The  validity  of  such  salinity  inferences  depends  completely  on  the  strength  of  the 
T-S  correlation,  or  what  Stommel  (1947)  called  "range  of  uncertainty",  i.c.,  the 
"tightness"  of  a  particular  T-S  cur\  e.  This  uncertainty  in  differential  salinity  can  be  re¬ 
lated  to: 

(17) 

where  AS  is  determined  from  the  T-S  curve,  r„  is  the  coefficient  of  determination  for 
temperature  and  salinity,  and  (1  —  rl)  is  the  fraction  of  salinity  variance  that  cannot  be 
explained  by  temperature  variance.  Clearly  the  uncertainty  increases  as  the  strength  of 
the  correlation  (and  magnitude  of »;,)  decreases.  Therefore,  it  may  be  reasonable  to  ex¬ 
pect  the  method  of  salinity  inference  from  T-S  relationship  to  be  most  accurate  in  re¬ 
gions  of  highest  T-S  correlations,  and  least  accurate  in  regions  where  the  T-S  correlation 
is  weak.  Blumberg  et  al.  (1984)  concluded  that  development  of  T-S  curves  for  the  CCS 
would  not  be  cost  clTectivc  for  modeling  efforts  due  to  combinations  of  relatively  poor 
T-S  correlations  inherent  in  near  surface  waters  (less  than  200  m),  and  lack  of  significant 
numbers  of  temperature  observations  in  deeper  areas  where  the  idea  might  otherwise 
work. 

Preliminary  investigations  have  been  made  to  determine  seasonal  values  of  the  T-S 
correlation  (i.e.,  coefficients  of  determination)  at  10,  100,  200  and  500  m.  Even  at 
500  m.  areas  where  temperature  variance  has  as  much  as  a  50  percent  linear  relationship 
with  salinity  are  rare  (Figure  62),  and  those  areas  are  not  constant  from  month  to 
month.  This  result  strongly  supports  tlie  statements  of  Blumberg  et  al.  (19S4).  A 
modeler  might  be  served  as  well  by  using  salinity  climatology  as  by  attempting  to  pre¬ 
scribe  tl'.c  salinity  field  based  on  its  relationship  with  temperature.  An 
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observationalist  should  probably  take  full  CTD  casts  to  determine  the  salinity  in  the 
CCS.  Calculations  could  be  performed  to  confirm  these  statements. 

The  geographic  distribution  of  salinity  in  the  CCS  may  be  sufficiently  linear,  at  least 
regionally,  to  allow  worthwhile  use  of  linear  regression  coefficients  on  latitude,  longitude 
and/or  depth.  Sums  required  by  least  squares  techniques,  to  compute  statistically  sig¬ 
nificant  coefficients  of  such  a  regression  equation  are  easily  obtainable  through  tech¬ 
niques  used  in  this  study.  Model  runs  using  the  resulting  salinity  parameterization 
would  be  instructive,  as  would  the  shorthand  description  of  the  salinity  field  represented 
by  the  parameterization  itself.  If  results  of  any  of  these  model  runs  differ  significantly 
from  corresponding  constant  salinity  runs,  decisions  concerning  the  best  method  to 
include  salinity  could  be  made. 

D.  SOUND  SPEED  CALCULATIONS 

Salinity  is  a  variable  in  the  equation  for  sound  speed  in  seawater  but  is  generally 
estimated  by  a  climatological  mean  when  tactical  acoustic  range  predictions  are  com¬ 
puted.  The  considerable  database  now  available  in  the  CCS  should  allow  a  numerical 
evaluation  of  the  validity  of  such  a  parameterization. 


91 


Figure  62.  Coefficient  of  Deterniinalion  (T-S),  Summer,  501)  m:  rj,  (the  algebraic 
sign  is  obtained  from  covariance,  and  "99"  indicates  a  bad  data  proc¬ 
essing  point). 


APPENDIX  A.  NUMBER  OF  OBSERVATIONS 

This  appendix  provides  plots  of  the  number  of  observations  associated  with  each 
data  block  (for  those  data  blocks  having  at  least  ten  observations)  for  the  months 
Januarj',  April,  July,  and  October  and  the  depths  10  and  500  m. 
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Figure  64.  Number  of  Observations  at  10  m,  April. 
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Figure  65.  Nuuil)cr  of  Oliscrvsilions  .il  10  lu,  .luly. 
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Figure  66.  Number  of  Observallons  at  10  m,  October. 
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Figure  67.  Number  of  Observ.itions  .il  500  m,  January. 


Figui'P  68.  Ntiinlier  of  Ohscrratioii.s  n(  50(1  in,  April. 
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Figure  70.  Number  «f  Ohscrvatious  at  500  lu.  October. 


APPENDIX  B.  SMALL  SCALE  EFFECTS  OF  SALINITY  ON  DYNAMIC 

HEIGHT 


In  addition  to  introducing  the  general  large-scale  biases  to  the  estimated 
geostrophic  flow  field  discussed  in  section  IVE,  the  estimator  may  misrepresent 
some  smaller  scale  features.  The  December  mean  map  shows  a  strong  closed  cyclonic 
circulation  near  30°  N,  121°  W.  The  December  mean  tbjsw  ®lso  shows  a  strong 
cyclonic  feature  at  that  point,  indicating  a  thermal  front  as  the  probable  cause  of  the 
enhanced  flow.  However,  the  difference  map  for  the  period  also  places  a  cyclonic  cir¬ 
culation  at  that  location,  showing  the  front  to  be  haline  as  well  as  thermal. 

In  April  a  closed  cyclone  again  appears  on  the  O'  map,  this  time  near  32°  N, 
123°  W.  Here  the  Oj.,,  map  shows  only  weak  cyclonic  turning  in  basically  onshore  flow. 
The  corresponding  O  map  indicates  strong  horizontal  cjxlonic  shear  in  the  region.  In 
this  case,  changes  in  salinity  over  a  relatively  small  distance  induce  considerable  shear 
in  the  mean  field.  This  feature  would  be  missed  entirely  by  an  estimate  based  on  tem¬ 
perature  variation  alone. 

In  May  a  similar  cyclonic  anomaly  appears  on  the  O'  map  near  36°  N,  127°  W. 
The  O  field  shows  a  strong  ridge  at  37°  30'  X  and  a  strong  trough  at  36°  X',  while  the 
Ojj,;.  map  shows  a  trough  at  37.5°X  and  a  ridge  at  36.5°X.  Admittedly  the  scale  of  this 
mean  feature  is  near  the  resolvable  limit,  and  some  critical  mean  values  of  O  have  been 
determined  from  relatively  few  casts.  However,  if  the  feature  is  correctly  represented,  it 
demonstrates  the  importance  of  salinity  as  a  modifier  of  the  effects  of  temperature.  It 
is  interesting  that  O  (200/500)  and  Oj,,^  (200/500)  contours,  by  contrast,  arc  in  phase  in 
this  area.  Analysis  of  the  May  salinity  distribution  over  the  area  of  interest  shows  sharp 
salinity  gradients  at  the  surface  and  less  pronounced  gradients  at  depth.  Consequently, 
the  inferred  flow  fields  are  consistent. 
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